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Abstract 

We present the current status of ongoing efforts to use functional methods, Dyson- 
Schwinger equations and functional renormalization group equations, for the descrip- 
tion of the infrared regime of nonabelian (pure) gauge theories in the Coulomb gauge. 
In particular, we present a new determination of the color-Coulomb potential with the 
help of the functional renormalization group that results in an almost linearly rising 
potential between static color charges at large spatial distances. 
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1 Introduction 



Important progress has been achieved over the last decade in the description of the deep 
infrared region of nonabehan gauge theories with the help of functional methods, employing 
Coulomb gauge fixing. By functional methods we refer to semi- analytical tools that do 
not make use of the discretization of space-time as does lattice gauge theory. Specifically, 
equations of Dyson-Schwinger type arising from a variational principle have been used, and 
more recently functional renormalization group equations. In this contribution, we will 
report on the current status of these investigations. We will focus exclusively on pure gauge 
theories, more specifically SU(N) Yang-Mills theory, but include static color charges so as to 
obtain a description of the heavy quark potential, as in quenched approximations. 

We will start by briefly describing the general theoretical setup: the Hamiltonian frame- 
work is used, where the Weyl and Coulomb gauge conditions, ^o(x) = and V • A"(x) = 0, 
are imposed on the SU(N) gauge fields. Physical states are described by wave functionals of 
A"(x) with scalar product 

= J D[A]J[A]<I>*[A]^[A]. (1) 

Here, J[A] stands for the Faddeev-Popov (FP) determinant J[A] = Dct (—V ■ D) with 
the spatial covariant derivative D"^ = 5°''V -|- gf/"^'^A'^(x). The functional integral in (1) is 
understood to be restricted to spatially transverse gauge fields, i.e., to those that fulfill the 
gauge fixing conditions. 

The dynamics is defined by the Christ-Lee Hamiltonian H [1] that we do not write out. 
In the presence of a static color charge density Pg(x), H contains the interaction term 

H, = IJ d'xd'yp-^{^)F^\^,y)pl{y) (2) 

with the integral kernel 

F"''(x, y) = (x, a|(-V • ■ D)-^|y, b) . (3) 

The vacuum expectation value (F"''(x, y)) is called the color-Coulomb potential. It is supp- 
posed to give the dominant contribution to the confining interaction between color charges. 
More precisely, for large spatial distances the color-Coulomb potential provides an upper 
bound for the Wilson potential [2]. 

For the following, it will be convenient to write the FP determinant in a local form by 
introducing ghost fields. 



J[A] = Det (-V • D) = y" D[c, c] exp j d^x c"(x)(-V • D"'')c^(x)^ . 



(4) 



In our analysis, we will focus on the equal-time correlation functions, i.e. the vacuum 
expectation values of products of the field operators A'*(x) (transverse), c"(x) and c"(x). 
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We can easily write down an expression for the generating functional of these correlation 
functions, 



X exp (^J (fx [J"(x) • A"(x) + c"(x)r;'^(x) + r;''(x)c"(x)]^ , (5) 

where V'fA] is the vacuum wave functional. If we formally define an "action" S[A\ through 
|V'[A]|^ = e"'^'''^', (5) looks like the usual generating functional of Euchdean Green's func- 
tions in the covariant Lagrangian formulation of the theory, only in three instead of four 
dimensions. Of course, S'[A] is a complicated and a priori unknown functional of A"(x). We 
will parametrize the "propagators", the equal-time two-point correlation functions of the 
theory, in the most general way (restricted by symmetries) as follows: 

(A^P) 4(-q)> = ^ (^^. - ^) (2^)''^(p - q) ' (6) 

(c"(p) c^(-q)> = <(p, a|(-V • D)-i|q, 6)> = ^ 5"" (27r)=^5(p - q) . (7) 

Here and in the following, we use the notation p = |p|. The functions uj{p) and d{p) will be 
of central interest in the rest of this contribution. Notice that the ghost propagator (7) is 
just the vacuum expectation value of the inverse FP operator (or rather, its integral kernel). 



2 Variational principle: Dyson- Schwinger equations 

A set of equations of Dyson- Schwinger type for the equal-time correlation functions of the 
theory was obtained in Ref. [3] from the variational principle, using a Gaussian ansatz for 
the vacuum wave functional. The contribution of the FP determinant was fully taken into 
account in [4,5]. We write the ansatz for the vacuum functional as 

|V^[A]|^ = e-^l^l, S[A] = \j ^^At{-v)2u{p)A1{^). (8) 

Then the variational principle with respect to the unknown function a;(p), 
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5u{p) 



(^) = , (9) 



leads to a gap equation for the equal-time gluon propagator. The detailed form of the 
equation as well as the approximations involved in its derivation can be found in [4,5]. 

The gap equation involves, apart from the gluon propagator, the ghost propagator and 
the color-Coulomb potential, hence further input is needed in order to arrive at a closed 
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system of equations. The generating functional (5) can be used to derive a Dyson- Schwinger 
(DS) equation for the ghost propagator in the usual way: 



p'd-\p) ^ • - = - ....iC!^-.-^ ■ (10) 



In the diagrams, we represent the full equal-time ghost propagator by a dashed line and 
the gluon propagator by a curly line, with a dot on the lines. By extending Taylor's non- 
renormalization theorem [6] to the present situation, we have replaced in (10) the full ghost- 
gluon vertex (one of the vertices on the right-hand side) with the bare one. This replacement 
is also used in the gap equation. 
For the color-Coulomb potential, 

{F^\p, -q)> = ((p, a|(-V • D)-i(-V2)(-V • D)-i|q, b)) = K(p) S^' (27r)=^5(p - q) , (11) 

we use the following parameterization and diagrammatic representation motivated by the 
appearance of the inverse FP operator [cf. (7)]: 

KM^f*pV(rtf^^ — , . (12) 

thereby defining the Coulomb form factor f{p). Of course, the function f{p) is itself un- 
known, and before discussing the possibility of determining it in terms of the gluon and ghost 
propagators in the next section, we will resort to the factorization hypothesis [7] 

((-V ■ D)-i(-V^)(-V • D)-i> = ((-V • D)-^> (-V^) ((-V • D)-i> , (13) 

which is equivalent to 

VM = ^P^'-^, (14) 

or f{p) = 1. Adopting this (so far unjustified) assumption, we obtain a closed system of 
equations. 

Before discussing the numerical solutions of the equations, we have to comment on the 
Gribov-Zwanziger confinement scenario [8,9]. In brief, the idea is that the existence of 
Gribov copies (gauge-equivalent but different transverse gauge field configurations) forces 
one to restrict the functional integral over the gauge field to the first Gribov region where the 
FP operator is positive definite.^ By a statistical argument, in the infrared (IR) regime the 
dominant contribution to the functional integral comes from the region close to the Gribov 



^Actually, the first Gribov region contains Gribov copies itself, and it is necessary to further restrict the 
integral to the so-called fundamental modular region. It is likely that this additional complication does not 
affect the following argument (in the case of the Coulomb gauge). 
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horizon where the FP operator has zero modes. Since the ghost propagator is the vacuum 

expectation value of the inverse FP operator, it may be argued that the ghost propagator 
d{p)/p^ should be more singular in the IR than p^"^, thus d~'^{p = 0) = 0, the "horizon 
condition" [10]. Hence, one should look for solutions that fulfill the horizon condition. 

It turns out that there are two different solutions of this type [4, 11].^ Both show scahng 
behavior in the IR, i.e., the equal-time propagators obey power laws in this kinematical 
regime. Furthermore, among the different contributions to the gluon propagator in the gap 
equation, the ghost loop diagram completely dominates the IR behavior, a property known 
as ghost dominance. These facts make it possible to even obtain analytical solutions for the 
propagators in the IR region [7, 12]. With the notations 

uj{p)^A p-" , d{p) = Bp-^ , (15) 

one obtains in this way a general sum rule for the IR exponents: 

a = 2/3 - 1 . (16) 

Consistent solutions exist for the values {a — 0.592,^ = 0.796) and (a = 1,^ = 1). One 
may also define a running coupling constant from the ghost-gluon vertex as 

"W-S^^. (17) 

[qb is the bare coupling constant). In the ultraviolet (UV), the solutions show asymptotic 
freedom (although not with the correct power of Inp due to the approximations to the gap 
equation), while a{p) saturates at a constant value in the IR. Analytically, one obtains 
Nca{Q) = 11.99 for the solution with /3 = 0.796, and N^aiO) = Wtt/S for /3 = 1 [12]. 

To close this section, we comment on a possible drastic simplification of the equations: if 
one uses, instead of the general Gaussian ansatz (8), the lowest-order perturbative vacuum 
wave functional 

iV'iAll^^e-^ol^l =exp (^-i J ^^^(-p)2pA"(p)) (18) 

in the generating functional (5), the complicated gap equation may be replaced by the 
following DS equation for the gluon propagator: 



-1 

2uj{p) = ^ = 2Z^p - -mrf frnir . (19) 



^Actually, in these numerical solutions a more sophisticated approximation of the Coulomb form factor, 
going beyond the factorization hypothesis, was used (see below). 
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Figure 1: Comparison of the results for the propagators from the two different sets of equa- 
tions [13]. The label "DSE" stands for the set that contains the gap equation and "FRG" 
for the set that uses (19). 

In particular, (19) does not make use of the factorization hypothesis. Due to ghost domi- 
nance, the gap equation approaches (19) in the JR. Somewhat surprisingly, we have found nu- 
merically that the solutions obtained with the two different sets of equations coincide over the 
whole momentum range from the IR to the UV to good numerical precision (for /3 = 0.796) 
[13]. In Fig. 1, ghost and gluon propagators are represented in a double-logarithmic plot 
for the two different sets of equations. The slight discrepancy in the IR is due to the lower 
numerical precision of the earlier calculation ("DSE") in [4]. 

3 Color- Coulomb potential and factorization hypothe- 
sis 

The color-Coulomb potential is more directly related to physically observable quantities than 
the gluon and ghost propagators. For the solution with /3 = 1, the factorization hypothesis 
(14) leads to Vc{p) oc j9~2-2/3 _ ^-4 ^j^g jj^ which corresponds to a potential in position 
space that rises exactly linearly for large distances. Unfortunately, the approximation used 
in the UV [see our remark following (17)] does not permit to relate the (Coulomb) string 
tension to the scale Aqcd- 

We will now turn to the question of whether the factorization hypothesis is actually 
justified. To this end, it is convenient to represent the color-Coulomb potential with the 
help of a composite operator K, 



((x,a|(-V-D)-i(-V2)(-V-D)-i|y 



b)) = (c«(x)irc^(y)> 



GI ' 




(20) 
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where the index GI (gluon-irreducible) on the vacuum expectation value means that one 

has to restrict the contributing diagrams to those where the operator K remains connected 
to the external points when all gluon lines are cut. The Coulomb form factor f{p) is then 
precisely the form factor of the composite operator K. Introducing K in the standard way 
in the generating functional (5), one may derive a DS equation for /(p). After suitable 
approximations, one arrives at (see also [4]) 



Now (21) can be used to close the system of equations instead of invoking the factorization 
hypothesis. The result is disappointing: no solution that fulfills the horizon condition could 
be found neither numerically nor analytically [14]. Numerically, solutions of the complete 
system of equations are found to exist only for (i~^(0) > 0.02. For the latter solutions, f{p) 
tends to a constant for p — )■ 0, so that Vc{p) oc and the color-Coulomb potential is not 
confining. We remark that this negative result is not due to the appearance of /(p) in the 
gap equation: in the simplified version that replaces the gap equation with (19), the two DS 
equations (10) and (19) decouple from (21), but still no solution can be found that would 
satisfy the horizon condition. We can now also explain the approximation that was used to 
obtain the numerical solutions [4, 11] (instead of the factorization hypothesis): the Coulomb 
form factor was determined from (21) with the ghost propagators in the loop diagram on 
the right-hand side replaced by the tree-level propagators. The result is a form factor that 
tends toward a constant in the JR. 

Let us briefly comment on the most recent results for the equal-time two-point correlation 
functions obtained in calculations on space-time lattices in the Coulomb gauge. They are 
still somewhat controversial. For the ghost propagator, a UV-bchavior d{p) oc (Inp)^^-^^^-'^) 
was found in [15] which is in reasonable agreement with the (also recent) result of one- 
loop perturbation theory d{p) oc (lnp)~°-^^ [16]. On the other hand, the IR-behavior was 
determined in [15] to be d{p) oc p-0-435(6) ^j-^j^h does not correspond to any of the two 
solutions of the DS equations. In [17], an IR-behavior uj{p) oc p^^ of the gluon propagator 
was found, which is consistent with the (/3 = l)-solution. For the UV-behavior of the 
gluon propagator, however, contradictory results have been reported: [17] finds u{p) oc p, as 
opposed to u{p) oc in [15]. Both results are at odds with the prediction of one- loop 

perturbation theory, uj{j)) oc p (lnp)°-^^ [16]. As for the color-Coulomb potential, although the 
recent results for the UV-behavior do not compare favorably with perturbation theory and 
the IR-behavior is inconclusive, it seems clear that the factorization hypothesis is violated 




(21) 



p 



[18]. 
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4 The functional renormalization group 

We will now turn to a different functional method, the functional (or Wilsonian) renormal- 
ization group. In order to adapt it to the case at hand, one starts with the generating 
functional (5) and introduces an IR cutoff k in the following way [13]: 



Z,[3,rj,fj]^ J D[c,c,A]exp(^- J ^ c«(-p)it:^(p) c"(p) 



The cutoff functions Rf.{p) and R^{p) have the properties 



RT^'P) ^ ^ for p<^k, Rk^ip) ^0 for p > 



(23) 



This means that the IR modes p ^ A: in the functional integral (22) arc heavily suppressed, 
while in the limit k ^ 0, R'jl^{p) — )■ and Zk[J,ri,f]] tends toward the full generating 
functional Z[J,rj,fj]. In the actual calculations, we have used an exponential suppression of 
the IR modes. 



RUp) = P'^rkip) , Rt{p) = 2prk{p) , rk{p) = exp 



P 



p_ 

,2 ^2 



(24) 



From (22), flow equations for the /c-dependent equal-time correlation functions can be derived 
in the standard way [19]. They read for the propagators 



2dkUk{p) = dk 



^,y'-Rtip) 



P 



(25) 



p'^dkd^^{p) = dk 



RUp) 




+ 



-♦-♦-<8>-»-5-<- 



(26) 



Here, the symbol 8) stands for the insertion of dkRk^- The non-renormalization theorem 

for the ghost-gluon vertex has been used in both equations. Furthermore, we have neglected 
diagrams that involve three- and four-gluon vertices, which is justified for the description of 
the IR regime if ghost dominance is assumed. Finally, we have omitted tadpole diagrams in 
order to be able to close the system of differential equations. Partial inclusion of the tadpole 
diagrams is argued in [13] to lead back to the DS equations (10) and (19) (after integrating 
over k). 
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Figure 2: The results for d{p) (left) and u!{p) (right) from the flow equations, for three 
different values of /cmin [13]. 

The general strategy is to start integrating the flow equations at a large value of k where 
due to asymptotic freedom the "action" S[A] can be replaced with 5*0 [A] from (18) and the 
coupling constant is small, so that the initial values of the flow are known. Then the flow 
equations are numerically integrated toward k = 0, where a;(p) = u!k=o{p) and d{p) = dk=Q{p) 
are read off. Technically, it is important to convert the differential equations (25), (26) to 
integral equations first, so that the horizon condition and a normalization condition for uj{p) 
can be conveniently incorporated. The results are presented in Fig. 2, again as double- 
logarithmic plots. For technical reasons, the integration of the flow equations stops at a 
minimum value kmm > 0. Then a;(p) = i^k=k^in{p) ior p ^ hnin, and similarly for d{p). From 
Fig. 2 it is clear that the power-law behavior of the propagators extends toward smaller 
momenta p as kmin is lowered. 

The exponents found numerically are (a = 0.28, f3 = 0.64), smaller than for both solutions 
of the DS equations. They obey the sum rule (16). The fact that the exponents come out 
smaller than the ones from the DS equations is not entirely unexpected, since a similar 
behavior was found in analogous calculations in the Landau gauge [20]. Generally, the 
results for the exponents will slightly vary with the choice of the cutoff functions due to the 
approximations made in the system of flow equations. An "optimized" choice is expected to 
give exponents identical to the ones from the DS equations [13]. For the running coupling 
constant (17) one also finds saturation in the IR at a slightly smaller value than for the DS 
solutions. 

By incorporating the composite operator K in the functional integral (22), one derives 
(after suitable approximations) the following flow equation for the Coulomb form factor: 

P'dkfk{p)=dj- ) =-^/^^^\^^ - - . (27) 

^ ' p p p 
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Figure 3: The Coulomb form factor f[p) from (27). 



Making use of the results for uJk(p) and dk{p), this equation can be integrated. Contrary 
to the DS equations, (27) has a solution that is represented in Fig. 3. The IR-behavior is 
determined numerically to 

fip)^p-\ 7 = 0.57. (28) 

In particular, of course, f{p) ^ 1, and the factorization hypothesis is violated. With the 
values for the exponents from the flow equations one obtains in the IR 

close to Vc{p) oc p~^ which would correspond to a linearly rising potential in position space. 

In summary, we find that functional methods are a powerful tool for the description of 
the nonperturbative infrared regime of nonabelian gauge theories. The formulation of these 
theories in the Coulomb gauge is particularly convenient, mainly because it gives direct ac- 
cess to the color-Coulomb potential. The Gribov-Zwanziger confinement scenario provides 
a conceptual framework to understand the confinement mechanism. It can be conveniently 
implemented via the horizon condition. In particular, we have seen that an almost lin- 
early rising color-Coulomb potential is obtained from the functional renormalization group 
equations (and the factorization hypothesis is violated). It has also become clear that the ap- 
proximations employed still have to be improved in order to achieve a quantitatively reliable 
description of the infrared region. 
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